CMBR Weak Lensing and HI 21-cm Cross-correlation Angular 

Power Spectrum 

Tapomoy Guha SarkaiQ 
Centre for Theoretical Studies, Indian Institute Of Technology. Kharagpur, 721302, India 

Abstract 

Weak gravitational lensing of the CMBR manifests as a secondary anisotropy in the temperature 
maps. The effect, quantified through the shear and convergence fields imprint the underlying large 
scale structure (LSS), geometry and evolution history of the Universe. It is hence perceived to be 
an important observational probe of cosmology. De-lensing the CMBR temperature maps is also 
crucial for detecting the gravitational wave generated B-modes. Future observations of redshifted 
21-cm radiation from the cosmological neutral hydrogen (HI) distribution hold the potential of 
probing the LSS over a large redshift range. We have investigated the correlation between post- 
reionization HI signal and weak lensing convergence field. Assuming that the HI follows the dark 
matter distribution, the cross-correlation angular power spectrum at a multipole i is found to be 
proportional to the cold dark matter power spectrum evaluated at i/r, where r denotes the comov- 
ing distance to the redshift where the HI is located. The amplitude of the cross-correlation depends 
on quantities specific to the HI distribution, growth of perturbations and also the underlying cos- 
mological model. In an ideal situation, we found that a statistically significant detection of the 
cross-correlation signal is possible. If detected, the cross-correlation signal holds the possibility of a 
joint estimation of cosmological parameters and also may be used to test various CMBR de-lensing 
estimators. 
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I. INTRODUCTION 

Weak gravitational lensing of distant background sources by intervening large scale 
structure, distorts their images over large angular scales. The effect arises due to fluctuations 
of the gravitational potential, and a consequent deflection of light by gravity. Measurement 
and quantitative study of these distortions allows us to probe the matter distribution and 
geometry of the universe. Late time evolution of the universe is dictated by dark energy 
through a modification of the growing mode of perturbations or through possible clustering 
properties of dark energy {w ^ —1). Weak lensing studies can be used to impose constraints 
on various cosmological parameters and hence, implicitly probe dark energy models ^ 
and modified gravity theories jsl. It is relevant for our present purpose to note that weak 
lensing is directly related to the underlying matter distribution of the universe. Weak 
lensing of background source galaxies by large scale structure (cosmic shear) has been studied 
extensively, and the measurements have been used for projected mass reconstruction (for 
review j4|) . 

Gravitational lensing also defiects the photons which are free streaming from the last scat- 
tering surface (epoch of recombination z ~ 1000) and manifests as a secondary anisotropy 
in the Cosmic microwave background radiation (henceforth CMBR) brightness tempera- 
ture maps 5[. Despite, the intrinsic weakness of the 'signal to noise ratio' for the above 
effect, weak lensing of CMBR can, in principle be used to probe the universe at distances 
{z ~ 1100) much larger than any galaxy- redshift surveys. Moreover CMBR lensing studies 
do not face the problems arising due to intrinsic alignment of source galaxies. Standard 
techniques to measure secondary anisotropies in CMBR, uses the cross correlation of some 
relevant observable (related to the CMB fluctuations) with fluctuations of some tracer of the 
large scale structure ^8\. Observables relevant to weak lensing are 'Convergence' and the 
'Shear' flelds, which quantify the distortion of an image due to gravitational lensing. Con- 
vergence (k) measures the lensing effect through its direct dependence on the gravitational 
potential and it probes geometry implicitly through its dependence on various cosmological 
distances. 

Future experiments (PLANK fl, CMBPOL js] etc) would provide high resolution maps 
for the CMB temperature and polarization flelds. The effect of gravitational lensing can 
be extracted from these maps by constructing various estimators for the convergence fleld 
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(k) through quadratic combination of these fields (T, E, B) |10l-|l2|. One could also predict 
the noise involved in such estimation based upon various experimental parameters. Lens- 
ing reconstruction can also be done using the 21 cm observations ISj. The reconstructed 
convergence field can then be used for cross correlation. De-lensing the CMB maps is also 
crucially important for detecting the gravitational wave generated B-mode. 

It is well accepted that the the neutral hydrogen (henceforth HI) distribution in the 
post-reionization epoch [z < 6) largely traces the underlying large scale structure of the 



universe 
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16| . This allows us to relate HI distribution to the cold dark matter distribution 



through a possible 'bias'. Matter perturbations are in the linear regime on large scales under 
consideration and the above simplifying assumption is reasonable. Hence, observations of 
the redshifted 21 cm radiation of the HI spin-fiip hyperfine transition provides an unique 



opportunity for probing the universe over a wide range of redshifts (200 > -2 > 0) 



14- 



16|. Theoretical predictions [l7|, [l8| have suggested the use of HI, statistically, as a probe 



of large scale _structure. Positive correlation between the optical galaxies (6dFGS) and HI 
fiuctuations [19-] has also been observed recently. 

In this paper we have investigated the possibility of using diffused cosmological HI as 
a tracer of the underlying large scale structure to probe weak lensing induced secondary 
anisotropy of the CMBR. Cosmic shear fields imprint the underlying distribution of matter 
over large scales. We have studied the the cross correlation between the post-reionization 
fiuctuations in the HI brightness temperature and the weak lensing convergence field. The 
cross-correlation angular power spectrum, measures the strength of the correlation as a 
function of the angular scale. 

The weak lensing of CMBR, quantified through the convergence field is expressed as a line 
of sight integral. Cross correlation of weak lensing with the HI fiuctuations, however pick 
up the contribution from only one redshift {zhi at which the HI is probed). The advantage 
of using HI observations is that, the redshifted 21 cm line emission observations allow us 
to probe the universe continuously at different redshifts. We can probe the integral effect 
of weak lensing at any intermediate redshift by suitably tuning the frequency band for HI 
observation. This, in principle enables us to do a tomographic study of the late-time cosmic 
history continuously over an entire range of redshifts. On similar lines, cross-correlation of 
HI temperature naap with the CMBR, aimed to isolate the ISW signal (an integral effect) 



has been studied 



20|. 
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Several Radio telescopes (eg. currently functioning GMRT [50| and upcoming MWA [51 1 



& LOFAR 



52| ) are aimed to map the cosmological distribution of HI at high redshifts. 



The extreme weakness of the post-reionization HI signal (< 10 /iJy) from individual clouds, 
despite some magnification due to Gravitational lensing [21'], poses a serious observational 
challenge. However, observation of the statistical distribution of HI as a weak background 
in radio observations does not require the need to resolve individual galaxies. Such observa- 
tions contain information about the HI fluctuations at the comoving distance being probed 
(frequency) 
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18|. 



Convergence field reconstructed from CMBR maps of large portion of the sky and a 
corresponding HI map would allow us to compute the cross-correlation power spectrum and 
hence independently quantify the cosmic history at redshifts z < 6. The cross-correlation 
power spectrum may also independently compare the various theoretical estimators that 
separate the lensing contribution from the CMB data . 



II. FORMULATION 

The lensed CMB brightness temperature T{n) along the direction of the unit vector n 
is related to the unlensed temperature T(n) through the map T(n) = T(n -|- a), where a 
denotes the total deflection due to weak lensing by the intervening large scale structure. At 
the lowest order, magnification of the signal is given by the convergence, k = — |V ■ a. The 
convergence field can be written as a line of sight integral given by [l| 

Q / H \ ^ rvLss 
^(n) = -n^o[-fj J drjFirj)S{dA{v)f^,v) (1) 

where (Ia stands for the comoving angular diameter distance and F{ri) is given by 

dAKriLSSjaym 

Here denotes the growing mode for the density contrast 5, and rj^ss denotes the confor- 
mal time corresponding to the last scattering surface (assuming instantaneous recombina- 
tion), and a(?7) denotes the scale factor. 

Here we have excluded weaker contribution to the convergence field from sources other 
than large scale structure (like gravitational waves). Expanding this in the basis of spherical 
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harmonics 

oo 

K{n) = J2(^emYUi^) (3) 

£,m 

The expansion coefficients a^^ can be obtained by integrating over the sohd angle as 



= / duJnK{n)Y;^{n) (4) 
Using the Raleigh expansion 

e^k - = 4nY,HYji{kr)Y;jmm{n) (5) 

£,m 

we have 

(-^)7 /o d^^(^)WjKM>^;™(k) (6) 

where 5(k) is the Fourier transform of S{r), and is the spherical Bessel function. 
In studying the post-reionization HI power spectrum we assume that the HI traces the 

1 /2 

underlying dark matter distribution with a possible bias function b{k) = [P^\k)/P{k)] , 
where P^\k) and P{k) denote the HI and dark matter power spectra respectively. This 
function is assumed to quantify the clustering property of the neutral gas. It is believed that, 
on small scales (below the Jean's length), the linear density contrast for the gas is related to 
the dark matter density contrast though a scale dependent function [22]. However the bias 
is known to be reasonably scale-independent on large scales. The length scale above which 
the bias is linear, depends crucially on the redshift being probed. Numerical simulations 



indicate that the large scale linear bias grows monotonically with redshift for 1 < 2; < 4 |23| 



This is known to be true for galaxies [24|-|26| . The increase in the amplitude of HI brightness 
temperature power spectrum is however slow (a factor of ~ 2 for z between 1 and 5) |27| . In 
this paper we have considered scales which are much larger than the scale of non-linearity 
and hence linear scale independent bias has been used. 



Expanding the HI 21-cm brightness temperature fluctuations (in Fourier space 28|]) from 
redshift zm in terms of spherical harmonics and proceeding as before we get 

——^5{\^^a)Ukr)YUk). (7) 



where xhi is the mean HI fraction, and 

0.02 J \ h J H{z) 



The term = k ■ n has its origin in the HI pecuhar velocities 17|, |28| which have also been 
assumed to be caused by the dark matter fluctuations. In equation ([7]) we have deflned 



Mx)=bje{x) - f 



(9) 



Where / denotes the logarithmic derivative of the growing mode and is given by / = fi! 



0.6 



At redshifts < 2 < 3.5 we have figas ~ 10 (for details see 29l-l31|). This allows us 



to calculate the mean neutral fraction of the hydrogen gas xhi = 50 f2gas^^(0.02/nb/;,^) = 
2.45 X 10^^, which we assume is a constant over the entire redshift range < 2; < 6. 

We use equations ([6]) and ([7]) to calculate Cf^~'^, the cross correlation angular power 
spectrum between the HI 21-cm brightness temperature signal and the convergence field, 
defined through 

(10) 



Hl-ni 



Note that Cf^~'^ also depends on zhi, the redshift from which the HI signal originates, or 
equivalently onv = 1420MHz/(l + 2;j|//), the frequency of the HI observations (not explicitly 
mentioned here). 
We obtain 



dk k'P{k)J,{krHi) / dvF{7^)je{kr 
where P{k) is the present day dark matter power spectrum, 

A{z) = ^Qmo(^] f{z)xHiD+{z) 



(11) 



(12) 

For large £ (small angular scales where "fiat sky" approximation is reasonable) the Limber 
approximation in Fourier space |32|, |33|, j^ikr) ^ ^/j^^oi^ + I ~ ^^)^ allows us to under- 
stand various generic scaling properties of the angular cross- correlation power spectrum. 



Cf, oc -A[zhi)— 

^ dAKZHl) 



(13) 



where P{k) is the present day dark matter power spectrum and all the terms on the rhs. 
are evaluated at zhi- 

Using equation ([6]) we have the Convergence auto-correlation power spectrum which for 
large ^ can be approximately written as 



dr]^A^P 



(14) 



We also have, for comparison, the HI-HI angular power spectrum Cf^zni) [35!], which 
describes the statistical properties of HI fluctuations . 

The function Cf^{zHi) is known to be a direct observational estimator of the HI fluctu- 
ations at redshift Zhi and does not require the assumption of an underlying cosmological 
model (eg. 3J]). Using the 'flat sky' approximation 35|, which is reasonable for I > 10, we 
have CI^\zhi) given by 

Cf\zHi) = —xl,Dl / dfc|| [b + f^^'] ' P{k) (15) 

'^'^u Jo 

where r is the comoving distance corresponding to the redshift z^i or equivalently frequency 
z/ = 1420MHz/(l + zm), and k = ^/c^ + (//r)2. In this paper we have used the WMAP5 
data for the various cosmological parameters. 

We note that the quantity of interest - the convergence field ^(n), is not a direct ob- 
servable in CMBR experiments. The degree of non-gaussianity in the lensed CMB maps is 
proportional to the lensing potential responsible for it. This allows a reconstruction of the 
weak lensing potential and consequently the deflection angle a, through the use of various 
statistical estimators {sg], [s^]. The reconstructed lensing convergence field is sensitive to the 
statistical tool (estimator) being used and reflects the degree of de-lensing achieved. 

The estimated quantity, we have focussed on, namely the cross correlation angular power 
spectrum, C^^"*^, does not directly de-lens the CMB maps. It however uses the recon- 
structed convergence field, and is hence sensitive to the underlying de-lensing technique, 
and the cosmological model. We have calculated the theoretical cross-correlation power 
spectrum assuming a standard cosmological model. The estimated Cf^~^, (where k is the 
estimated convergence field) with its known error bars can be compared with our predicted 
Cf^~^. Hence, the theoretical cross-correlation angular power spectrum provides a template 
to independently compare various estimators which are aimed at de-lensing the CMB maps. 



III. RESULTS 



Figure [T] shows the theoretically predicted cross-correlation angular power spectrum 
qHI-k various redshifts 0.5 < zhi < 6. The currently favored ACDM cosmological 
model with parameters (fimo, ^ao, ^5 cts, '^s) = (0.28, 0.72, 0.7, 0.82, 0.97) [38|, 



391] has been 



used here. For HI signal we have assumed a linear bias model (reasonable on the large scales 
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FIG. 1: The HI- convergence angular power spectrum for redshifts z = 0.5.1.0,2.0 and 5.0 (top to 
bottom). 



under consideration) with 6 = 1 in the fiducial model. 

Numerical simulations have indicated the deviation from 6 = 1 at high redshifts. It is seen 
that at large scales the linear bias is 6 ~ 2 for 2 ~ 3. The effect of larger (scale independent) 
bias is shown in figure |2l Apart from the scaling of the power spectrum at large scales 
the bias also has a weak effect of modifying the power spectrum amplitude through the 
change in the the redshift space distortion factor f3 = f /b. We have also indicated the scale 
^NL ~ ^NL^^z, above which the linear bias assumtion is invalid. For 2; ~ 3 this angular scale 
/nl ~ 6000. We have restricted ourselves to multipoles less than that. 

Figure |3] shows the Convergence auto-correlation power spectrum for reference. The 
Cross-correlation power spectrum has the same shape as the matter power spectrum. For 
different redshifts the signal peaks at a particular i which scales as £ oc vhi- The angular 
distribution of power clearly follows the underlying clustering properties of matter. The 
amplitude of the cross-correlation power spectrum depends on various factors some of which 
are related to the underlying cosmological model and others related to the HI distribution 
at zhi- The angular diameter distances directly imprint the geometry of the universe and 
also depends on the cosmological parameters. The 21cm signal has been proposed to be 
an useful probe of the cosmological parameters 
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42| . The cross-correlation signal may 
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FIG. 2: The HI angular power spectrum at redshift z = 3.3 showing the effect of hnear bias, /nl 
is the scale above which non-linear biasing should be incorporated 



likewise be used independently for joint estimation of parameters. 

We shall now discuss the prospect of detecting the cross-correlation signal. Redshifted 
21 cm signal is buried deep under foregrounds. Removal of the foreground component is a 



major challenge 
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43(1 . However, it is to be noted that cross-correlation between the 



HI brightness temperature field and the convergence field is much less likely to be affected 
by foregrounds or other systematics. The error in the cross-correlation signal is a sum in 
quadrature, of the contribution due to instrumental noise and sample variance. Increased 
resolution (for CMB experiment) and increased time of observation (for 21 cm observation) 
can in principle significantly reduce the instrumental noise. Sample variance however puts 
a fundamental bound on the detectability of the signal. 

The sample variance for the cross- correlation angular power spectrum Cf^~^ is given by 

4 = ^li}^ (16) 

Where the numerator contains the auto-correlation angular power spectra. A£ represents 
a band in i and fs is fraction of sky common to the convergence field k and HI bright- 
ness temperature distribution denotes the number of independent estimates of the 
21cm observations obtained from different frequency channels in a given frequency band and 
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FIG. 3: The convergence power spectrum 



suppresses the sample variance by a factor 1 / \/Wc- 

We have used the ideal hypothetical possibility of = 1, and used A£ = 1 . we have 
chosen = 32 assuming that the HI signal decorrelates over a frequency separation of ~ 
IMHz and hence yield 32 independent estimates for a 32MHz bandwidth radio observation. 
The estimated Signal to Noise ratio S/N = Cf^~'^/asv is shown in figlH . for zhi = 0.5. 
The predicted S/N is seen to be ~ 2 and is not high enough for a statistically significant 
detection which requires S/N > 3. Choosing a Ai = 10 for i < 100 and Ai = 100 for 
i > 100 will however produce a S/N > 3. 

It is possible to increase the S/N by collapsing the signal from different scales i and 
thereby test the feasibility of a statistically significant detection. The Signal to Noise cumu- 
lated upto a multipole i is defined as (see \^ for similar calculation) 



^(2£+l)v^/,(C, 



HI-K\ 



{cr+Nrm+N^ ^^^^ 

The summation in the above equation extends up to a certain i. N^ and Nf^ denotes the 
noise power spectrum for k and HI observations respectively. Ignoring the instrument noises 
we note that there is a significant increase in the S/N by cumulating over multipoles £. 
This implies that a statistically significant detection of Cf^~^ is possible and the signal is 
not cosmic variance limited. 21- cm observations allow us to probe a continuous range of 
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FIG. 4: The lower curve shows Signal to Noise ratio (S/N) as a function of angular scale I. The 
upper curve shows the effect of summing over multipoles. The probing redshift zui = 0.5 



redshifts. This allows us to further increase the SjN by collapsing the signal from various 
redshifts. As discussed earlier, an increased HI bias would increase the signal. However the 
SjN oc C^^^'^ j \fC^ is not expected to be seriously affected. 

Instrumental noise plays an important role at large multipoles (small scale). For a typ- 



ical CMB experiment, the noise power spectrum 45|, |46| is given by = ap-^Jl^pi^Wi ^, 



where different pixels are assumed to have uncorrelated noise with uniform variance 



a 



pix 



s^/^pix, where and tpix denotes pixel sensitivity and 'time spent on the pixel' 
respectively, ilpix is the solid angle subtended per pixel and we choose a gaussian beam 
We = exp[— £^^^|^jjp^/161n2]. For CMBPOL [9| like experiments, the noise power spectrum 
for K with the beam FWHM ~ 3' and sensitivity ~ lui^ — arcmin is < 10~^ for i < 3000 

nn 

(see [9|, |47|) Hence, <^ and maybe ignored in our present analysis. 

For HI observations, the quantity of interest is the complex Visibility which is used to 
estimate the power spectrum 34j. For a radio telescope with N antennae, system temper- 
ature Tsys, operating at a frequency z/, and band width B the noise correlation is given by 



48| 



HI 



sys 



1 



2ttN{N-1) [ K J TVA^UO ^'Am ^'piU^iy) 



(18) 
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Where 27rf/ ~ ^, T denotes total observation time, and K is related to the effective 
collecting area of the antenna dish . The function p{U, v) takes any non-uniform distribution 
of baselines into account and depends on the array design. The bin AU = I/ttOq is chosen 
assuming a gaussian beam of width 6q. With a GMRT or MWA like instrument js^, one 
can in principle achieve a noise level much lesser than the signal by increasing the time of 
observation (a 2000 hour observation is sufficient even with the present GMRT cofiguration) 
and also by increasing the band width of the instrument. Being inversely related to the 
number of antennae in the array, future designs can allow further suppression of the the 
system noise and achieve N^^ « Cf^. 

This establishes the detectability of the cross-correlation signal. We would like to conclude 
by noting that this theoretical prediction of positive correlation between weak lensing fields 
and 21 cm maps, quantified through Cf^~'^ may allow an independent means to estimate 
various cosmological parameters and also test various estimators for CMBR delensing. 
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CMBR Weak Lensing and HI 21-cm Cross-correlation Angular 

Power Spectrum 

Tapomoy Guha SarkaiQ 
Centre for Theoretical Studies, Indian Institute Of Technology. Kharagpur, 721302, India 

Abstract 

Weak gravitational lensing of the CMBR manifests as a secondary anisotropy in the temperature 
maps. The effect, quantified through the shear and convergence fields imprint the underlying large 
scale structure (LSS), geometry and evolution history of the Universe. It is hence perceived to be 
an important observational probe of cosmology. De-lensing the CMBR temperature maps is also 
crucial for detecting the gravitational wave generated B-modes. Future observations of redshifted 
21-cm radiation from the cosmological neutral hydrogen (HI) distribution hold the potential of 
probing the LSS over a large redshift range. We have investigated the correlation between post- 
reionization HI signal and weak lensing convergence field. Assuming that the HI follows the dark 
matter distribution , the cross-correlation angular power spectrum at a multipole I is found to 
be proportional to the cold dark matter power spectrum evaluated at i/r, where r denotes the 
comoving distance to the redshift where the HI is located. The amplitude of the cross-correlation 
depends on quantities specific to the HI distribution, growth of perturbations and also the underly- 
ing cosmological model. In an ideal situation, we found that a statistically significant detection of 
the cross-correlation signal is possible. If detected, the cross-correlation signal hold the possibility 
of a joint estimation of cosmological parameters and also test various CMBR de-lensing estimators. 

Keywords: gravitational lensing; intergalactic media; power spectrum 
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I. INTRODUCTION 

Weak gravitational lensing of distant background sources by intervening large scale 
structure distorts their images over large angular scales. The effect arises due to fluctuations 
of the gravitational potential, and a consequent deflection of light by gravity. Measurement 
and quantitative study of these distortions allow us to probe the matter distribution and 
geometry of the universe. Late time evolution of the universe is dictated by dark energy 
through a modification of the growing mode of perturbations or through possible clustering 
properties of dark energy {w ^ —1). Weak lensing studies can be used to impose constraints 
on various cosmological parameters and hence, implicitly probe dark energy models ^ 
and modified gravity theories jsl. It is relevant for our present purpose to note that weak 
lensing is directly related to the underlying matter distribution of the universe. Weak 
lensing of background source galaxies by large scale structure (cosmic shear) has been studied 
extensively, and the measurements have been used for projected mass reconstruction (for 
review j4|) . 

Gravitational lensing also defiects the photons which are free streaming from the last scat- 
tering surface (epoch of recombination z ~ 1000) and manifests as a secondary anisotropy 
in the Cosmic microwave background radiation (henceforth CMBR) brightness tempera- 
ture maps 5[. Despite, the intrinsic weakness of the 'signal to noise ratio' for the above 
effect, weak lensing of CMBR can, in principle be used to probe the universe at distances 
{z ~ 1100) much larger than any galaxy- redshift surveys. Moreover CMBR lensing studies 
do not face the problems arising due to intrinsic alignment of source galaxies. Standard 
technique to measure secondary anisotropics in CMBR, uses the cross correlation of some 
relevant observable (related to the CMB fluctuations) with fluctuations of some tracer of 
large scale structure Is-Sj. Observables relevant to weak lensing are 'Convergence' and the 
'Shear' flelds, which quantifles the distortion of an image due to gravitational lensing. Con- 
vergence (k) measures the lensing effect through its direct dependence on the gravitational 
potential and it probes geometry through its dependence of various cosmological distances. 

Future experiments (PLANK H, CMBPOL jg] etc) would provide high resolution maps 
for the CMB temperature and polarization flelds. The effect of gravitational lensing can 
be extracted from these maps by constructing various estimators for the convergence fleld 
(k) through quadratic combination of these flelds (T, E, B) |10l-ll2|. One could also predict 



2 



the noise involved in such estimation based upon various experimental parameters. Lens- 



ing reconstruction can also be done using the 21 cm observations [13j. The reconstructed 
convergence field can then be used for cross correlation. De-lensing the CMB maps is also 
crucially important for detecting the gravitational wave generated B-mode. 

It is well accepted that the the neutral hydrogen (henceforth HI) distribution in the 
post-reionization epoch {z < 6) largely traces the underlying large scale structure of the 



universe 



1J-[16| . This allows us to relate HI distribution to the cold dark matter distribution 



through a possible 'bias'. Matter perturbations are in the linear regime on large scales under 
consideration and the above simplifying assumption is reasonable. Hence, observations of 
the redshifted 21 cm radiation of the HI spin-flip hyperfine transition provides an unique 



opp ortunity for probing the universe over a wide range of redshifts (200 > z > 0) 
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16|. Theoretical predictions [l7|, [l8| have suggested the use of HI, statistically, as a probe 



of large scale structure. Positive correlation between the optical galaxies (6dFGS) and HI 
fiuct nations has also been observed recently. 

In this paper we have investigated the possibility of using diffused cosmological HI as 
a tracer of the underlying large scale structure to probe weak lensing induced secondary 
anisotropy of the CMBR. Cosmic shear fields imprint the underlying distribution of matter 
over large scales. We have studied the the cross correlation between the post-reionization 
fiuctuations in the HI brightness temperature and the weak lensing convergence field. The 
cross-correlation angular power spectrum, measures the strength of the correlation as a 
function of angular scale. 

The weak lensing of CMBR, quantified through the convergence field is expressed as 
a line of sight integral. Cross correlating this with HI fiuctuation however picks up the 
contribution from only one redshift [zhi at which the HI is probed). The advantage of 
using HI observations is that, the redshifted 21 cm line emission observations allow us to 
probe the universe continuously at different redshifts. We can probe the integral effect of 
weak lensing at any intermediate redshift by suitably tuning the frequency band for HI 
observation. This, in principle enables us to do a tomographic study of the late-time cosmic 
history continuously over an entire range of redshifts. On similar lines, cross-correlation of 
HI temperature rnap with the CMBR, aimed to isolate the ISW signal (an integral effect) 



has been studied 



2Q|. 



Several Radio telescopes (eg. currently functioning GMRT 39| and upcoming MWA 40 1 



3 



& LOFAR 



41| ) are aimed to map the cosmological distribution of HI at high redshifts. 



The extreme weakness of the post-reionization HI signal (< 10 /iJy) from individual clouds, 
despite some magnification due to Gravitational lensing [21], poses a serious observational 
challenge. However, observation of the statistical distribution of HI as a weak background 
in radio observations however does not require the need to resolve individual galaxies. Such 



observation contain in: 
probed (frequency) 117 



brmation about the HI fluctuations at the comoving distance being 



18|- 



Convergence field reconstructed from CMBR maps of large portion of the sky and a 
corresponding HI map would allow us to compute the cross-correlation power spectrum and 
hence independently quantify the cosmic history at redshifts z < 6. The cross-correlation 
power spectrum may also independently compare the various theoretical estimators that 
separate the lensing contribution from the CMB data . 



II. FORMULATION 

The lensed CMB brightness temperature T(n) along the direction of the unit vector fi 

is related to the unlensed temperature T(n) through the map T(n) = T(n + a), where oc 

denotes the total deflection due to weak lensing by intervening large scale structure. At 

the lowest order, magnification of the signal is given by the convergence, k, = — |V ■ a. The 

convergence field can be written as a line of sight integral given by [l| 

Q / H \ ^ rvLss 
/^(n) = -^]„o(^^J J dr]F{r])6{dA{r])h,r]) (1) 

where cIa stands for the comoving angular diameter distance and F{ri) is given by 

^ dAiVLSs - v)dA{r])D+{v) 
dA{VLSs)a{'n) 

Here denotes the growing mode for the density contrast 6, and rjiss denotes the confor- 
mal time corresponding to the last scattering surface (assuming instantaneous recombina- 
tion), and a{ri) denotes the scale factor. 

Here we have excluded weaker contribution to the convergence field from sources other 
than large scale structure (like gravitational waves). Expanding this in the basis of spherical 
harmonics 

fi:(n) = ^<„F,„(n) (3) 

£,m 



The expansion coefficients can be obtained by integrating over the sohd angle ua as 

«L = j duf,Kin)Y;^in) (4) 

Using the Raleigh expansion 

47r5^(-0^j,(fcr)y4(k)r,„(n) (5) 



we have 

aL = 6vrr]^o(^^ J H)'y dry F(r/)5(k)j,(A:r)r;jk) (6) 

where 6(h) is the Fourier transform of 5(r), and je{x) is the spherical Bessel function. 
In studying the HI powers spectrum we assume that the HI traces the underlying dark 
matter distribution with a possible bias b. 



Expanding the HI 21-cm brightness temperature fluctuations (in Fourier space 22|) from 
redshift zhi in terms of spherical harmonics and proceeding as before we get 

5{k,a)Mkr)Y;M- (7) 

(27r) 



where xhi is the mean HI fraction, and 

0.02 J \ h J H{z) 



f(.)=4.0mK(l + .)^(5^U^)^ (8) 



The term = k • n has its origin in the HI peculiar velocities 17|, |22| which have also been 



assumed to be caused by the dark matter fluctuations. In equation ([7]) we have defined 

J,(x)=bUx) - . (9) 



At redshifts < z < 3.5 we have Qg^s ~ 10~'^ (for details see 23N25|). This allows us 
to calculate the mean neutral fraction of the hydrogen gas xhi = 50 f2gas/i^(0.02/i7{,/i^) = 
2.45 X 10~^, which we assume is a constant over the entire redshift range < 2; < 6. 

We use equations (jS]) and ([7]) to calculate Cf^~^, the cross correlation angular power 
spectrum between the HI 21-cm brightness temperature signal and the convergence field, 
defined through 

('^£m'^£'m') = '^^ee'Smm' (10) 

Note that Cf^~'^ also depends on zhi the redshift from which the HI signal originates, or 
equivalently on u = 1420MHz/(l + zhi) the frequency of the HI observations, but we do 
not show this explicitly here. 



We obtain 



C 



Hl-n 



--A{zhi) 



dk 



eP{k)Je{krHi) 



VLSS 



driF{ri)je{kr) 



vo 



(11) 



where P{k) is the present day dark matter power spectrum, 



TT 



A{z) = -^mo{ — T{z)xhiD+{z 



(12) 



For large i (small angular scales where "flat s^ 



ty" approximation is reasonable) we can use 



the Limber approximation in Fourier space 26|, |27|], which allows us to replace the spherical 
Bessel functions je{x) by a Dirac deltas Sd{x). We have ji{kr) ^ y^j^^oi^-^ + I ~ ^'^)' 
whereby the angular cross-correlation power spectrum takes the simple form 



^lA[zni)ih + f)I^P 



(13) 



2 ' ' •" dA{zHl) 

where P{k) is the present day dark matter power spectrum and all the terms on the rhs. 
are evaluated at zhi. 

Using equation ([6]) we have the Convergence auto-correlation power spectrum which for 
large ^ can be approximately written as 

We also have, for comparison, the HI-HI angular power spectrum Cf^{zHi) [29^, which 
describes the statistical properties of HI fluctuations . 

The function C^\zhi) is known to be a direct observational estimator of the HI fluctu- 
ations at redshift Zm and does not require the assumption of an underlying cosmological 
model (eg. 28|). Using the 'flat sky' approximation 29|, which is reasonable for i > 10, we 
have Cj^\zHi) given by 

f2 



Cf\zHl) 



—xl,Dl dk\\[b + ff,'YP{k) 



(15) 



where r is the comoving distance corresponding to the redshift zm or equivalently frequency 



u = 1420MHz/ (1 + zm), and k = ^kj + {l/rf. 

In this paper we have used the WMAP5 data for the various cosmological parameters. 
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FIG. 1: The HI- convergence angular power spectrum for redshifts z = 0.5.1.0,2.0 and 5.0 (top to 
bottom). 
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FIG. 2: The convergence power spectrum C? 



III. RESULTS 



Figure [T] shows the theoretically predicted cross-correlation angular power spectrum 
qHI-k various redshifts 0.5 < zhi < 6. The currently favored ACDM cosmological 
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model with parameters {Qrno,^Ao,h,as,ns) = (0.28, 0.72, 0.7, 0.82, 0.97) [30|, l3l| has been 
used here. For HI signal we have assumed a linear bias model (reasonable on the large 
scales under consideration) with 6 = 1 in the fiducial model. Figure [2] shows the Conver- 
gence auto-correlation power spectrum for reference. The Cross-correlation power spectrum 
has the same shape as the matter power spectrum. For different redshifts the signal peaks 
at a particular i which scales as £ oc thi- The angular distribution of power clearly fol- 
lows the underlying clustering properties of matter. The amplitude of the cross-correlation 
power spectrum depends on various factors some of which are related to the underlying 
cosmological model and others related to the HI distribution at zhi- The angular diameter 
distances directly imprint the geometry of the universe and also depends on the cosmological 
parameters. The 21cm signal has been proposed to be an useful probe of the cosmological 



parameters 32H34j|. The cross-correlation signal may likewise be used independently for 
joint estimation of parameters. 

We shall now discuss the prospect of detecting the cross-correlation signal. Redshifted 



21 cm signal is buried 



deep under foregrounds. Removal of the foreground component is a 



major challenge 28|, l29|, |34, |35|| . However, it is to be noted that cross-correlation between the 
HI brightness temperature field and the convergence field is much less likely to be affected 
by foregrounds or other systematics. The error in the cross-correlation signal is a sum in 
quadrature, of the contribution due to instrumental noise and sample variance. Increased 
resolution (for CMB experiment) and increased time of observation (for 21 cm observation) 
can in principle significantly reduce the instrumental noise. Sample variance however puts 
a fundamental bound on the detectability of the signal. 



The sample variance for the cross-correlation angular power spectrum Q is given by 



{2(i + l)^JM 

Where the numerator contains the auto-correlation angular power spectra. A£ represents 
a band in ^ and fs is fraction of sky common to the convergence field k and HI bright- 
ness temperature distribution Nf. denotes the number of independent estimates of the 
21cm observations obtained from different frequency channels in a given frequency band and 
suppresses the sample variance by a factor 1 / ^/Wl. 

We have used the ideal hypothetical possibility of = 1, and used A£ = 1 . we have 
chosen iVc = 32 assuming that the HI signal decorrelates over a frequency separation of ~ 
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100 




1 

FIG. 3: The lower curve shows Signal to Noise ratio (S/N) as a function of angular scale I. The 
upper curve shows the effect of summing over multipoles. The probing redshift zui = 0.5 



IMHz and hence yield 32 independent estimates for a 32MHz bandwidth radio observation. 
The estimated Signal to Noise ratio SjN = Cf^~'^/asv is shown in fig[3] . for zhi = 0.5. 
The predicted S/N is seen to be ~ 2 and is not high enough for a statistically significant 
detection which requires S/N > 3. Choosing a Ai = 10 for i < 100 and A£ = 100 for 
i > 100 will however produce a S/N > 3. 

It is possible to increase the S/N by collapsing the signal from different scales i and 
thereby test the feasibility of a statistically significant detection. The Signal to Noise cumu- 
lated upto a multipole i is defined as (see [36|] for similar calculation) 



S_yj^i2i+l)^JsiCf'- 



(17) 

N^ and N!^^ denotes 



Nj ^(ci^i+Nrm+m 

The summation in the above equation extends up to a certain 
the noise power spectrum for n and HI observations respectively. Ignoring the instrument 
noises we note that there is a significant increase in the S/N by cumulating over multipoles 
£. This implies that a statistically significant detection of Cf^~'^ is possible and the signal 
is not cosmic variance limited. 21cm observations allow us to probe a continuous range of 
redshifts. This allows us to further increase the S/N by collapsing the signal from various 
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redshifts. 

Instrumental noise plays an important role at large multipoles (small scale). For 
CMBPOL {9] like experiments, the noise power spectrum for k for i < 1000 is roughly 
flat and considerably less than the signal (see jo], Is^]). For HI observations with a GMRT 



like instrument [28|] , one can in principle achieve a noise level much lesser than the signal by 
increasing the time of observation and also by increasing the band width of the instrument. 
Being inversely related to the number of antennae in the array, future designs can allow 
further suppression of the the system noise. 

This establishes the detectability of the cross-correlation signal. We would like to conclude 
by noting that this theoretical prediction of positive correlation between weak lensing fields 
and 21 cm maps, quantified through Cf^~'^ may allow an independent means to estimate 
various cosmological parameters and also test various estimators for CMBR delensing. 
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